Fluorescence-based directed termination PCR: direct mutation characterization without sequencing

We describe a fluorescence-based directed termination PCR (fluorescent DT–PCR) that allows accurate determination of actual sequence changes without dideoxy DNA sequencing. This is achieved using near infrared dye-labeled primers and performing two PCR reactions under low and unbalanced dNTP concentrations. Visualization of resulting termination fragments is accomplished with a dual dye Li-cor DNA sequencer. As each DT–PCR reaction generates two sets of terminating fragments, a pair of complementary reactions with limiting dATP and dCTP collectively provide information on the entire sequence of a target DNA, allowing an accurate determination of any base change. Blind analysis of 78 mutants of the supF reporter gene using fluorescent DT–PCR not only correctly determined the nature and position of all types of substitution mutations in the supF gene, but also allowed rapid scanning of the signature sequences among identical mutations. The method provides simplicity in the generation of terminating fragments and 100% accuracy in mutation characterization. Fluorescent DT–PCR was successfully used to generate a UV-induced spectrum of mutations in the supF gene following replication on a single plate of human DNA repair-deficient cells. We anticipate that the automated DT–PCR method will serve as a cost-effective alternative to dideoxy sequencing in studies involving large-scale analysis for nucleotide sequence changes.

Single-pair fluorescence resonance energy transfer on freely diffusing molecules: Observation of Förster distance dependence and subpopulations

Photon bursts from single diffusing donor-acceptor labeled macromolecules were used to measure intramolecular distances and identify subpopulations of freely diffusing macromolecules in a heterogeneous ensemble. By using DNA as a rigid spacer, a series of constructs with varying intramolecular donor-acceptor spacings were used to measure the mean and distribution width of fluorescence resonance energy transfer (FRET) efficiencies as a function of distance. The mean single-pair FRET efficiencies qualitatively follow the distance dependence predicted by Förster theory. Possible contributions to the widths of the FRET efficiency distributions are discussed, and potential applications in the study of biopolymer conformational dynamics are suggested. The ability to measure intramolecular (and intermolecular) distances for single molecules implies the ability to distinguish and monitor subpopulations of molecules in a mixture with different distances or conformational states. This is demonstrated by monitoring substrate and product subpopulations before and after a restriction endonuclease cleavage reaction. Distance measurements at single-molecule resolution also should facilitate the study of complex reactions such as biopolymer folding. To this end, the denaturation of a DNA hairpin was examined by using single-pair FRET.

Polymerase recognition of synthetic oligodeoxyribonucleotides incorporating degenerate pyrimidine and purine bases

A universal base that is capable of substituting for any of the four natural bases in DNA would be of great utility in both mutagenesis and recombinant DNA experiments. This paper describes the properties of oligonucleotides incorporating two degenerate bases, the pyrimidine base 6H,8H-3,4-dihydropyrimido[4,5-c][1,2]oxazin-7-one and the purine base N6-methoxy-2,6-diaminopurine, designated P and K, respectively. An equimolar mixture of the analogues P and K (called M) acts, in primers, as a universal base. The thermal stability of oligonucleotide duplexes were only slightly reduced when natural bases were replaced by P or K. Templates containing the modified bases were copied by Taq polymerase; P behaved as thymine in 60% of copying events and as cytosine in 40%, whereas K behaved as if it were guanine (13%) or adenine (87%). The dUTPase gene of Caenorhabditis elegans, which we have found to contain three nonidentical homologous repeats, was used as a model system to test the use of these bases in primers for DNA synthesis. A pair of oligodeoxyribonucleotides, each 20 residues long and containing an equimolar mixture of P and K at six positions, primed with high specificity both T7 DNA polymerase in sequencing reactions and Taq polymerase in PCRs; no nonspecific amplification was obtained on genomic DNA of C. elegans. Use of P and K can significantly reduce the complexity of degenerate oligonucleotide mixtures, and when used together, P and K can act as a universal base.

Ultraviolet plumage colors predict mate preferences in starlings

Avian plumage has long been used to test theories of sexual selection, with humans assessing the colors. However, many birds see in the ultraviolet (<400 nm), to which humans are blind. Consequently, it is important to know whether natural variation in UV reflectance from plumage functions in sexual signaling. We show that female starlings rank males differently when UV wavelengths are present or absent. Principal component analysis of ≈1300 reflectance spectra (300–700 nm) taken from sexually dimorphic plumage regions of males predicted preference under the UV+ treatment. Under UV− conditions, females ranked males in a different and nonrandom order, but plumage reflectance in the human visible spectrum did not predict choice. Natural variation in UV reflectance is thus important in avian mate assessment, and the prevailing light environment can have profound effects on observed mating preferences.

What initiates speciation in passion-vine butterflies?

Studies of the continuum between geographic races and species provide the clearest insights into the causes of speciation. Here we report on mate choice and hybrid viability experiments in a pair of warningly colored butterflies, Heliconius erato and Heliconius himera, that maintain their genetic integrity in the face of hybridization. Hybrid sterility and inviability have been unimportant in the early stages of speciation of these two Heliconius. We find no evidence of reduced fecundity, egg hatch, or larval survival nor increases in developmental time in three generations of hybrid crosses. Instead, speciation in this pair appears to have been catalyzed by the association of strong mating preferences with divergence in warning coloration and ecology. In mate choice experiments, matings between the two species are a tenth as likely as matings within species. F1 hybrids of both sexes mate frequently with both pure forms. However, male F1 progeny from crosses between H. himera mothers and H. erato fathers have somewhat reduced mating success. The strong barrier to gene flow provided by divergence in mate preference is probably enhanced by frequency-dependent predation against hybrids similar to the type known to occur across interracial hybrid zones of H. erato. In addition, the transition between this pair falls at the boundary between wet and dry forest, and rare hybrids may also be selected against because they are poorly adapted to either biotope. These results add to a growing body of evidence that challenge the importance of genomic incompatibilities in the earliest stages of speciation.

A pair of adjacent glucocorticoid response elements regulate expression of two mouse metallothionein genes

Synthesis of mouse metallothionein (MT)-I and MT-II is transcriptionally induced by the synthetic glucocorticoid, dexamethasone (DEX) or both in vivo as well as in numerous cell lines. However, the location(s) of a glucocorticoid response element (GRE) has not been described. The observation that a marked MT-I gene, as well as heterologous genes, when placed in the context of 17 kb of flanking sequence from the MT locus, are inducible by DEX and lipopolysaccharide in transgenic mice renewed the search for the GRE. Analysis of a series of deletion constructs from this 17-kb region in cultured cells identified a single 455-bp region that conferred DEX induction on a reporter gene. This 455-bp region contains two GREs that bind to the glucocorticoid receptor as assessed by gel mobility shift. Deletion of this fragment from the 17-kb flanking region eliminates the DEX responsiveness of reporter genes. The two GREs, which are located ≈1 kb upstream of the MT-II gene and ≈7 kb upstream of the MT-I gene, are necessary for induction of both genes and can function independently of elements within the proximal promoter region of either gene.

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

NMR scalar couplings across Watson–Crick base pair hydrogen bonds in DNA observed by transverse relaxation-optimized spectroscopy

This paper describes the NMR observation of 15N—15N and 1H—15N scalar couplings across the hydrogen bonds in Watson–Crick base pairs in a DNA duplex, hJNN and hJHN. These couplings represent new parameters of interest for both structural studies of DNA and theoretical investigations into the nature of the hydrogen bonds. Two dimensional [15N,1H]-transverse relaxation-optimized spectroscopy (TROSY) with a 15N-labeled 14-mer DNA duplex was used to measure hJNN, which is in the range 6–7 Hz, and the two-dimensional hJNN-correlation-[15N,1H]-TROSY experiment was used to correlate the chemical shifts of pairs of hydrogen bond-related 15N spins and to observe, for the first time, hJHN scalar couplings, with values in the range 2–3.6 Hz. TROSY-based studies of scalar couplings across hydrogen bonds should be applicable for large molecular sizes, including protein-bound nucleic acids.

A genome-wide search for susceptibility genes in human systemic lupus erythematosus sib-pair families

Systemic lupus erythematosus (SLE) is an autoimmune multisystem inflammatory disease characterized by the production of pathogenic autoantibodies. Previous genetic studies have suggested associations with HLA Class II alleles, complement gene deficiencies, and Fc receptor polymorphisms; however, it is likely that other genes contribute to SLE susceptibility and pathogenesis. Here, we report the results of a genome-wide microsatellite marker screen in 105 SLE sib-pair families. By using multipoint nonparametric methods, the strongest evidence for linkage was found near the HLA locus (6p11-p21) [D6S257, logarithm of odds (lod) = 3.90, P = 0.000011] and at three additional regions: 16q13 (D16S415, lod = 3.64, P = 0.000022), 14q21–23 (D14S276, lod = 2.81, P = 0.00016), and 20p12 (D20S186, lod = 2.62, P = 0.00025). Another nine regions (1p36, 1p13, 1q42, 2p15, 2q21–33, 3cent-q11, 4q28, 11p15, and 15q26) were identified with lod scores ≥1.00. These data support the hypothesis that multiple genes, including one in the HLA region, influence susceptibility to human SLE.

Evolution of paired domains: Isolation and sequencing of jellyfish and hydra Pax genes related to Pax-5 and Pax-6

Pax proteins are a family of transcription factors with a highly conserved paired domain; many members also contain a paired-type homeodomain and/or an octapeptide. Nine mammalian Pax genes are known and classified into four subgroups: Pax-1/9, Pax-2/5/8, Pax-3/7, and Pax-4/6. Most of these genes are involved in nervous system development. In particular, Pax-6 is a key regulator that controls eye development in vertebrates and Drosophila. Although the Pax-4/6 subgroup seems to be more closely related to Pax-2/5/8 than to Pax-3/7 or Pax-1/9, its evolutionary origin is unknown. We therefore searched for a Pax-6 homolog and related genes in Cnidaria, which is the lowest phylum of animals that possess a nervous system and eyes. A sea nettle (a jellyfish) genomic library was constructed and two pax genes (Pax-A and -B) were isolated and partially sequenced. Surprisingly, unlike most known Pax genes, the paired box in these two genes contains no intron. In addition, the complete cDNA sequences of hydra Pax-A and -B were obtained. Hydra Pax-B contains both the homeodomain and the octapeptide, whereas hydra Pax-A contains neither. DNA binding assays showed that sea nettle Pax-A and -B and hydra Pax-A paired domains bound to a Pax-5/6 site and a Pax-5 site, although hydra Pax-B paired domain bound neither. An alignment of all available paired domain sequences revealed two highly conserved regions, which cover the DNA binding contact positions. Phylogenetic analysis showed that Pax-A and especially Pax-B were more closely related to Pax-2/5/8 and Pax-4/6 than to Pax-1/9 or Pax-3/7 and that the Pax genes can be classified into two supergroups: Pax-A/Pax-B/Pax-2/5/8/4/6 and Pax-1/9/3/7. From this analysis and the gene structure, we propose that modern Pax-4/6 and Pax-2/5/8 genes evolved from an ancestral gene similar to cnidarian Pax-B, having both the homeodomain and the octapeptide.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

A novel pair of immunoglobulin-like receptors expressed by B cells and myeloid cells

An Fcα receptor probe of human origin was used to identify novel members of the Ig gene superfamily in mice. Paired Ig-like receptors, named PIR-A and PIR-B, are predicted from sequence analysis of the cDNAs isolated from a mouse splenic library. Both type I transmembrane proteins possess similar ectodomains with six Ig-like loops, but have different transmembrane and cytoplasmic regions. The predicted PIR-A protein has a short cytoplasmic tail and a charged Arg residue in the transmembrane region that, by analogy with the FcαR relative, suggests the potential for association with an additional transmembrane protein to form a signal transducing unit. In contrast, the PIR-B protein has an uncharged transmembrane region and a long cytoplasmic tail containing four potential immunoreceptor tyrosine-based inhibitory motifs. These features are shared by the related killer inhibitory receptors. PIR-A proteins appear to be highly variable, in that predicted peptide sequences differ for seven randomly selected PIR-A clones, whereas PIR-B cDNA clones are invariant. Southern blot analysis with PIR-B and PIR-A-specific probes suggests only one PIR-B gene and multiple PIR-A genes. The PIR-A and PIR-B genes are expressed in B lymphocytes and myeloid lineage cells, wherein both are expressed simultaneously. The characteristics of the highly-conserved PIR-A and PIR-B genes and their coordinate cellular expression suggest a potential regulatory role in humoral, inflammatory, and allergic responses.

Nucleic acid duplexes incorporating a dissociable covalent base pair

We have used molecular modeling techniques to design a dissociable covalently bonded base pair that can replace a Watson-Crick base pair in a nucleic acid with minimal distortion of the structure of the double helix. We introduced this base pair into a potential precursor of a nucleic acid double helix by chemical synthesis and have demonstrated efficient nonenzymatic template-directed ligation of the free hydroxyl groups of the base pair with appropriate short oligonucleotides. The nonenzymatic ligation reactions, which are characteristic of base paired nucleic acid structures, are abolished when the covalent base pair is reduced and becomes noncoplanar. This suggests that the covalent base pair linking the two strands in the duplex is compatible with a minimally distorted nucleic acid double-helical structure.

Proper sorting of the cation-dependent mannose 6-phosphate receptor in endosomes depends on a pair of aromatic amino acids in its cytoplasmic tail

The 67-amino acid cytoplasmic tail of the cation-dependent mannose 6-phosphate receptor (CD-MPR) contains a signal(s) that prevents the receptor from entering lysosomes where it would be degraded. To identify the key residues required for proper endosomal sorting, we analyzed the intracellular distribution of mutant forms of the receptor by Percoll density gradients. A receptor with a Trp19 → Ala substitution in the cytoplasmic tail was highly missorted to lysosomes whereas receptors with either Phe18 → Ala or Phe13 → Ala mutations were partially defective in avoiding transport to lysosomes. Analysis of double and triple mutants confirmed the key role of Trp19 for sorting of the CD-MPR in endosomes, with Phe18, Phe13, and several neighboring residues contributing to this function. The addition of the Phe18-Trp19 motif of the CD-MPR to the cytoplasmic tail of the lysosomal membrane protein Lamp1 was sufficient to partially impair its delivery to lysosomes. Replacing Phe18 and Trp19 with other aromatic amino acids did not impair endosomal sorting of the CD-MPR, indicating that two aromatic residues located at these positions are sufficient to prevent the receptor from trafficking to lysosomes. However, alterations in the spacing of the diaromatic amino acid sequence relative to the transmembrane domain resulted in receptor accumulation in lysosomes. These findings indicate that the endosomal sorting of the CD-MPR depends on the correct presentation of a diaromatic amino acid-containing motif in its cytoplasmic tail. Because a diaromatic amino acid sequence is also present in the cytoplasmic tail of other receptors known to be internalized from the plasma membrane, this feature may prove to be a general determinant for endosomal sorting.